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AN IMPROVED PROCESS FOR PHASE-DERIVED RANGE MEASUREMENTS 
BACKGROUND 

The invention relates generally to radar systems, and more particularly, to radar range 
5 measurement. 

Typically, a radar system (or more simply a radar) tracks a target's centroid position to 
obtain an estimate of ballistic flight characteristics for range determination. To obtain estimates 
of ballistic flight, radar systems typically use filters such as Kalman filters of relatively high 
order to produce range measurement data. Such filtering is designed to model only a target's 
10 ballistic trajectory. Therefore, it is possible to track and measure range from target trajectory 
quite accurately. 

To observe small variations in range due to motions other than ballistic motion, it is 
necessary to sample target returns in a manner that ensures that phase change measured from one 
radar pulse to the next is unambiguous. Such unambiguous phase measurements require a 
1 5 significant amount of radar system resources, more specifically, higher signal-to-noise (SNR) 

and data rates, thus reducing the safety margin built into the radar to cope with interferers such as 
main beam jamming or an increased number of targets. 

SUMMARY 

20 The present invention features a technique that provides a high precision range 

measurement at reduced signal-to-noise ratio (SNR) and data rate so that radar system objectives 
in a main beam jamming and/or multiple target environment can be achieved. 

In one aspect, therefore, determining range of a radar target includes receiving signal 
samples based on returns of a target during tracking, processing the signal samples to produce a 

25 wideband envelope range estimate for components of target motion comprising precession and 
spin components, measuring the signal samples to produce ambiguous phase values and 
using each wideband envelope range estimate and ambiguous phase value to produce an 
unambiguous phase value indicative of range. 

Embodiments of the invention may include one or more of the following features. 

30 The range determination can further include determining an estimate of ballistic 

trajectory of a target for the signal samples and removing the estimated ballistic trajectory during 
processing. 
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The processing of the signal samples can operate at a sampling rate that is at least twice 
the frequency of the target spin motion components. 

Using each wideband envelope range estimate and ambiguous phase value to produce the 
unambiguous phase value can include subtracting the measured ambiguous phase from the 
5 wideband envelope range estimate to produce an error value associated with the wideband 
envelope range estimate and subtracting the error value from the wideband envelope range 
estimate to give the unambiguous phase value. 

The magnitude of the error value can be determined and resources of a radar system that 
performs the tracking can be adjusted to ensure that the magnitude of the error is sufficiently 
10 low. 

The signal sample processing can include: i) producing a spectrum of wideband 
envelope range estimates from the signal samples; ii) transforming the wideband envelope range 
estimates to obtain a spectral estimate of each motion component of precession, spin, spin plus 
precession and spin minus precession; iii) detecting each motion component; iv) estimating 

15 amplitude, frequency and phase for each motion component spectral estimate; and v) forming a 
sinusoid in range motion from the estimate of amplitude, frequency and phase for each motion 
component spectral estimate. 

The processing can occur in batch mode for signal samples obtained during several 
cycles of precession motion. 

20 The sinusoid in range motion can be used to determine an integer number k of cycles in 

phase change between pulses of the signal samples. A value of 2rck can be added to the 
measured ambiguous phase value prior to subtracting the measured ambiguous phase value from 
the wideband envelope range estimate. 

In another aspect, a system includes a transmitter/receiver to direct transmit signals to and 

25 receive return signals from a target, a processor to process the return signals as in-phase and 
quadrature samples to produce angle information and range signals, a tracker to track a target 
detected according to results of the processing by the first processor, the tracker measuring range 
data during tracking and estimating a ballistic trajectory therefrom, and a unit operable to use a 
wideband envelope range estimate of the in-phase and quadrature samples, the ballistic trajectory 

30 estimate and an ambiguous phase measurement of the in-phase and quadrature samples to 
produce a range measurement that is unambiguous in phase. 
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Particular implementations of the invention may provide one or more of the following 
advantages. Through the use of the range phase determination of the present invention it is 
possible to observe small motions in body dynamics of targets, e.g. re-entry vehicles and decoys, 
and measure range of such targets with the precision of phase. Moreover, it is possible to reduce 
5 the data rate and signal-to-noise ratio requirements of phase-derived range measurement by using 
wideband envelope range (WBER). 

Other features and advantages of the invention will be apparent from the following 
detailed description and from the claims. 

1 o BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of a tracking radar system that includes a data processor which 
performs phase-derived range measurements. 

FIG. 2 is a block diagram of an exemplary data processor that employs a phase-derived 
range (PDR) unit to produce unambiguous phase values for phase-derived range measurements. 
1 5 FIG. 3 is a plot of range vs. time for an exemplary ballistic trajectory of a target as well 

as the associated precession and spin components of target motion. 

FIG. 4 is a block diagram of a PDR unit. 

FIG. 5 is a flow diagram illustrating the operation of a Wideband Envelope Range 
(WBER) estimator of a PDR unit. 
20 FIG. 6 is a plot of range versus time. 

FIG. 7 is a plot of phase versus range. 

FIG. 8 is a plot of spectrum of range relative to center of gravity. 
FIG. 9 shows a plot of spectral signal-to-noise requirements. 
FIG. 10 is a plot of precession range estimation errors. 
25 FIG. 1 1 is a plot of spin range estimation errors. 

FIG. 12 is a plot of precession phase error differences. 
FIG. 1 3 is a plot of spin phase error difference. 

FIG. 14 is a plot of phase difference of 1Hz Sine wave for several data rates. 
FIG. 15 is a plot illustrating the tradeoff between SNR and data rate for constant power. 
30 Like reference numerals will be used to represent like elements. 
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DETAILED DESCRIPTION 

Referring to FIG. 1, a tracking radar system 10 is shown. The tracking radar system 10 
may be a ground-based radar system, but could be employed on a ship or in an aircraft or other 
airborne or non-airborne vehicle as well. In one embodiment, the tracking radar system 10 is 

5 implemented as a pulse Doppler tracking radar system. The tracking radar system 10 includes a 
transmitter 12, the output of which is delivered to an antenna system 14 for radiation in the form 
of a transmit beam. The antenna system 14 collects echo (return) signals received from a target. 
The return signals are processed by a receiver 16 to detect the presence of the target and 
determine its location in range and in angle. The antenna system 14 can include a mechanically 

10 or electronically rotating antenna. A duplexer 1 8 coupled to the transmitter 12, receiver 16 and 
antenna system 14 allows the antenna system 14 to be used on a time-shared basis for both 
transmitting and receiving. 

Still referring to FIG. 1, the receiver 16 includes a low-noise amplifier ("LNA") 20 that is 
coupled to a down converter 22, which performs radio frequency (RF) to intermediate frequency 

15 (IF) conversion. A receiver exciter or local oscillator ("LO") 24 provides oscillator signals to the 
down converter 22. The down converter 22 is connected to a signal processor 26 via an AID unit 
28. The echo signals as received by the signal processor 26 correspond to in-phase ("I") and 
quadrature ("Q") signals. The signal processor 26 operates on the I and Q signals (as sampled 
by the A/D converter 28) to produce processed I, Q samples 30 and angle information 32, e.g., 

20 elevation and/or azimuth angle information. The signal processor 26 performs filtering, possibly 
including pulse compression filtering, envelope detection and post-detection integration (video 
integration), and angle determination, among other functions. 

The system 1 0 further includes a threshold detector 34 and a data processor 36. The 
threshold detector 34 receives the I, Q samples from the signal processor 26 and generates from 

25 those samples a range (gate) signal 38. It will be appreciated that the threshold detection and 
signal processing may be executed by a single processor or computer. The signal processor 26 
and threshold detector 34 can be coupled to, and provide detection information to other 
conventional radar system elements, e.g., the data processor 36, as well as a display and control 
center, e.g., a ballistic management control center (not shown). The I, Q samples produced by 

30 the signal processor 26 and range signal (from the threshold detector 34) are provided to the data 
processor 36 for use in further processing, as will be described. Data processor results (output 
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40) may be provided to the control center, as indicated. Those aspects of the system 10 and its 
components not described herein can be implemented according to known techniques. 

Referring to FIG. 2, the data processor 36 includes a tracker 50, a phase-derived range 
(PDR) block or unit 52 and a classification unit 54. The tracker 50 receives the gated I/Q 

5 samples 30 and angle information 32 from the signal processor 26, and receives the range gate 
signal 38 from the threshold detector 34. The tracker 50 initiates a track for a detected target, 
conducts tracking to measure range data and estimates from the range data a target trajectory 55 
(using a target trajectory estimator 56). The tracker 50 provides to the PDR block 52 range gated 
I/Q samples 58 for the tracked target. The PDR block 52 operates on these samples to produce a 

10 precession and spin spectral analysis, as will be described later with reference to FIG. 4. The 
results of this spectral analysis (body dynamics 60) are provided to the classification unit 54, 
which uses this and other information (e.g., output 62 from the tracker 50 and angle information 
32 from the signal processor 26) to provide a comprehensive set of target characteristics 64 to the 
control center (not shown). The characteristics can include, for example, information that 

15 indicates how a target is to be classified, for example, classification as threatening or non- 
threatening. The control center uses this information (and, possibly, information from other 
sources, such as satellites) to determine if the target classification is sufficient or if tracking 
should continue. The control center sends commands to the radar transmitter/receiver to indicate 
which of the targets are a priority, and thus types of waveforms, data rates and SNR to use. 

20 The tracker 50 initiates a target track, conducts tracking for a predetermined time period, 

e.g., between 25 -100 seconds, and stores range data measured on each pulse. The tracker 50 
batch-processes the entire batch of range data with a relatively high-order filter that accurately 
predicts ballistic trajectory. More specifically, the tracker 50 applies the high-order filter to the 
range data, in batch mode, to estimate a ballistic trajectory of the target. The target ballistic 

25 trajectory estimate is removed (during subsequent PDR block processing) to leave only the target 
precession and spin motion components, the phase of which can be estimated with sufficient 
precision, as will be described below with reference to FIG. 4. 

FIG. 3 illustrates, in a plot of range versus time 70, the various components of target 
motion. These motion components include target ballistic path or trajectory 72, spin 74 and 

30 precession 76. The spin and precession motions follow a sine wave type of pattern, as shown in 
the figure. Precession refers to wobbling (or lateral displacement) motion imparted to a spinning 
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target during flight. For the example illustrated, the target motion has a precession motion with 
a frequency of 1/8 Hz and peak amplitude of 0. 1 m. The spin motion for the illustrated example 
has a spin motion frequency of 1 Hz and a peak amplitude of .01 m. The minimum Nyquist 
sampling at which the precession motion can be sampled (assuming the ballistic trajectory 
5 component has been removed) is twice the precession frequency. However, it is the spin 
frequency that determines how often the precession motion must be sampled. Thus, the 
minimum sampling frequency is twice the spin motion frequency, or 2 Hz for the illustrated 
example. 

Referring to FIG. 4, the PDR block 52 receives from the tracker 50 the gated I, Q samples 

10 58 and the ballistic trajectory estimate 55. The PDR block includes a Wideband Envelope Range 
(WBER) estimator 80 and a unit 82 for measuring or obtaining measurements of phase "A" 84 of 
the I, Q samples. Both of these units receive and operate on the I, Q samples 58. The measured 
phase 84 is precise but ambiguous in that the number of cycles is unknown. The WBER 
estimator 80 produces from the samples a WBER phase estimate 86 in the form "27tk + A + e", 

15 where "k" is the integer number of cycles in phase change from one pulse to the next, and "e" is 
the WBER estimate error value. The error needs to be less than ± n radians for k to be an integer 
number. The unit 52 further includes a "k" estimator 88, and a "27rk" adder 90 to convert the k 
cycles to radians and add the "2n:k" radians to the ambiguous phase A to form a phase value 92 
of "27tk + A". The PDR unit 52 includes two subtractors, a first subtractor 94 and a second 

20 subtractor 96. The first subtractor 94 subtracts the phase value 92 from the WBER phase 

estimate 86 to produce an error value 98 indicative of the error e. The second subtractor 96 
subtracts the error value 98 (that is, "s") from the WBER phase estimate 86 to produce an 
unambiguous phase value 100 indicative of range. The resulting output 27tk+A, which gives the 
range with the precision of phase, is unambiguous as long as the value of the WBER error e is 

25 less than 1 radian root-mean-squared (RMS). It will be understood that one or more of the 
elements represented in the block diagram may be implemented in software. 

The SNR and length of dwell time are adjusted so that the error (in residue) pulse-to- 
pulse is small compared to ± n radians (or la, approximately 1 radian). A goal of lo error of 1 
radian is selected so that ±3a = +3 radians (~ n radians) ensures that range, with the precision of 

30 phase, can be measured unambiguously. 
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The approach and assumptions for the spectral analysis of the PDR unit 52 are as 
follows. The pulse repetition frequency (PRF) that is used is just above the Nyquist sampling 
rate, that is, e.g., equal to or greater than twice the spin frequency. The processing of the PDR 
unit 52 can be achieved by batch processing, rather than with a real-time phase-locked loop, so 
5 that the required data rate is determined by Nyquist sampling as opposed to phase-locked loop 
requirements. 

A high order filter in the WBER estimator 80 batch-processes several cycles of precession 
motion. The output of the filter contains a smooth estimate of range. Wideband pulses are used 
to estimate envelope range. A spectrum of wideband envelope range estimates produce 

10 estimates of precession, spin and spin ± precession motion for long dwells and high SNR. In the 
analysis, it is assumed that WBER estimates give the range to a target's center of gravity. That 
is, the smoothed range of the filter is varying slowly, so that pulse-to-pulse differences between 
range measured and estimated by the filter is small. 

Referring to FIG. 5, the WBER estimator 80 operates as follows. The WBER estimator 

15 80 receives as input the I, Q samples (step 90). The WBER estimator 80 produces a spectrum of 
WBER estimates from the input samples (step 92). The WBER estimator 80 removes the 
ballistic trajectory estimated by the tracker (step 94). The WBER estimator 80 Fourier- 
transforms the WBER estimates to obtain a spectral estimate of each motion component of 
precession, spin and spin ± precession (step 96). The resolved precession and spin components 

20 are detected (if their peak excursions are greater than XI A) (step 98), and their amplitude a p , 
frequency f, and phase 0 are measured (step 100). The WBER estimator 80 constructs a 
sinusoid in range r = ap sin (27ift +0) from the spectral estimates of amplitude ap, frequency f and 
phase 0 for each motion component (step 1 02). The WBER estimator 80 computes the 
differential dr to find amplitude, frequency and phase errors in range (step 104). Referring back 

25 to FIG. 4, the integer value of "k" is determined once the frequency estimate is known. 

FIGS. 6-15 illustrate results of WBER estimator processing for an example with the 
following numerical values: 

30 f s = 1 Hz spin frequency 

fp = 1/8 Hz precession frequency 
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Np = 3 precession cycles 
PRF = 3 Hz 

S/N = 15 dB, per-pulse SNR 

X = .03 m, wavelength 

Signal bandwidth = 2B 

Rp = . 1 m, peak precession range 

R s = .01 m, peak spin range 

R s .p = .0025 m, peak spin minus precession range 

R s +p = .007 m, peak spin plus precession range 

The precession and spin components have peak precession and peak spin values representative of 
actual targets. The spin-precession component is plotted with higher than expected value, to 
illustrate its position in the spectral plots. The spin plus precession component can be larger than 
the spin component, depending on trajectory and target motion scenarios. 

Referring to FIG. 6, a plot of range in the time domain for three precession cycles (24 
seconds) shows that precession dominates, with the spin and other components modulating at 1 - 
.125 = .875 Hz, 1 Hz, 1+.125 = 1.25 Hz. The data is sampled at 3 Hz, Fourier transformed and 
zero padded, and then inverse Fourier transformed to obtain a near continuous plot. 

Referring to FIG. 7, a plot of phase versus range is shown. This plot gives peak phase 
excursion in cycles for a given peak range in meters. For example, for a 0. 1 meter peak the 
number of cycles is: 

.1 m x (4n radians/^m) = 41.9 radians = 6.67 cycles. 
The spin component is .01 m peak or 4.19 radian = .67 cycles, or just slightly more than +.5 
cycle. Thus, the spin component in the example requires detection with sufficient spectral SNR 
to accurately model the amplitude, phase and frequency. As discussed earlier, components with 
a peak excursion of less than .0075 m need not be detected and modeled. 

Referring to FIG. 8, a plot of the spectrum of WBER samples (relative to the target 
centroid position) is shown. The plot is the discrete Fourier transform (DFT) of the WBER 
measurements for 3 cycles of precession (8 second precession period) and a 1 Hz spin frequency. 
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It will be appreciated that if only spin were present, 2 Hz (twice the spin frequency) 
constitutes Nyquist sampling. Such a sampling frequency is twice the usual radar frequency and 
is the consequence of the data being purely real. 

In the example, a 3 Hz sampling rate was selected to account for spin plus precession 

5 components and to provide sampling margin. A number of different techniques may be used to 
select sampling frequency. One option is to use a priori data. Another option would be to apply 
a bootstrap technique. 

The minimum data dwell is determined by spacing between spectral lines. Precession 
lines adjacent to the spin line are separated from the spin line by 1/T P = .125 Hz, the precession 

10 frequency. The dwell should have two or three periods to provide good resolution, particularly 
when weighting is applied. The plot shows that light weighting should be added to avoid 
detection of sidelobes in the vicinity of the spin line. The RMS noise line shows noise about -36 
dB relative to precession and about 13 dB below the spin line. Thus, the spin line can also be 
detected. The 3 Hz data rate therefore provides a spectrum in which all the useful data lies 

15 between 0 and 1.5 Hz. 

The precession line for the illustrated example is sufficiently large that it can be readily 
detected. For the spin components, the critical requirement of SNR, data rate and dwell time is 
to provide a spectral line large enough so that spin can be detected, and with sufficiently low 
error that the data is effectively unwrapped in phase. It is not necessary to achieve a detection of 

20 the spectral line for small components. As long as the phase excursion a p is sufficiently small so 

6,71X1 p <n or a. <- = .0075 meters, or less than + Vz cycle, the output remains phase 
X 4 

unambiguous. The peak excursion a p is in meters, and ~- times the peak excursion results in 

range phase in radians. 

The amplitude is estimated from the magnitude of the line. For each line, the frequency 
25 is estimated from linear odd weighting of the time data to obtain a "monopulse" estimate in the 
frequency domain. Phase is estimated from the ratio of quadratic to real component. Noise 
errors for each estimate of amplitude, frequency, and phase are uncorrelated with each other, 
although not independent. The differential dr gives the error s in each of the three components, 
amplitude, frequency and phase, that is, e a , £f and £e, respectively, as described earlier. 
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Amplitude, frequency and phase errors in range are Gaussian, zero mean, uncorrelated 
and with standard deviations: 

<r a = sm(2rft + 0)* Eq. 1 
V 25 / 

a f = 2nta p co&tft + 0) * A { Eq. 2 

a B = a p cos(2^ + 0) * -jL= Eq. 3 

V 2, V 

where 

\33cycles 

AK = 

IB 

Af = \/{N s -\)r 

N s = number of samples 



S f = N s \ ~y | , spectral signal-to-noise ratio Eq. 4 

B = wideband chirp bandwidth 
r = time between samples 

Results using Eqs. 1 -3 are consistent with the estimation of sinusoidal signals techniques 
as described in "Fundamentals of Statistical Signal Processing: Estimation Theory," by Steven 
M. Kay (University of Rhode Island, Prentice Hall, 1993). 

From the expression for spectral signal-to-noise ratio S f (Eq. 4) it can be seen that S f 
depends on peak amplitude of the sinusoidal motion and bandwidth. 

Referring to FIG. 9, which shows a plot illustrating spectral SNR (S f ) requirements, it is 
observed that, for a signal slightly less than Log (a p ) = -2, no spectral detection is required for 
this or smaller values. Thus, the requirement for SNR and signal bandwidth is to produce a 
spectral value of S f (slightly greater than 10 dB, depending on number of samples and allowable 
false alarm rate) adequate for the desired P d . 

Referring to FIG. 10, a plot of precession WBER estimation errors as a function of time 
is shown. The plot gives the la variation of each of the amplitude, frequency and phase errors 
separately. Except for the frequency error at the edge of the window, it should be noted that the 
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errors are less than +.5 radian. With WBER range error less than ±.5 radians, the number of 
cycles k can be estimated so that the error of ±.5 radians is sufficient to avoid ambiguities. Thus, 
the RMS value of the sum of the estimation values at any specific time is less than the 1 radian 
goal. Weighting of the data to reduce sidelobes can further de-emphasize the edge values. Near 
5 the center of the window, frequency errors are smaller than the amplitude and phase components. 

Referring to FIG. 1 1, a plot of spin range estimation error as a function of time is shown. 
The time reference in the plot is arbitrary. The plot shows the effect separately of amplitude, 
frequency, and phase errors on the reconstruction of the spin line. It can be noted that the 
amplitude and phase errors are the same magnitude, but shifted by 90°. This follows from 



FIG. 12 shows a plot of precession phase error differences as a function of time. The 
low frequency of the precession motion produces phase differences of negligible amplitude 
compared to +1 radian. 

FIG. 13 shows a plot of spin phase error difference as a function of time. At 3 Hz and 15 
15 dB S/N, the spin phase difference would be approximately 1 .7 times the results shown in FIG. 

11. For this reason, the plot of FIG. 14 shows the results of operating at 10 Hz with 10 dB SNR. 

Referring to FIG. 14, consider an error sine wave of amplitude 1 radian (as shown). If 
this sine wave is sampled at the Nyquist rate (2 Hz), the phase difference peaks at ±2 radians. 
For higher sampling rates the phase difference reduces, as shown in the figure. For example, at 3 
20 Hz, a 1 .7x increase over the 4 Hz example is observed. Note that at 6 Hz the phase difference is 
+1 radians. These results suggest that a tradeoff between SNR and data rate will result in a 
different balance, where higher data rate with low SNR is favored by a phase difference loop or a 
phase difference Doppler, compared to a batch process where data rate much higher than Nyquist 
is not needed. 

25 FIG. 1 5 shows a plot of data rate (PRF) versus SNR, and illustrates the tradeoff between 

these two resources for constant power. The curve is a plot of SNR x data rate = a constant. For 
example, 15 dB SNR at a PRF of 3 Hz would require the same energy as a 12 dB SNR at a PRF 
of 6 Hz. For a batch process with a 3 Hz PRF, a 15 dB SNR is preferable for higher P d in the 
spectral detection process. 
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Eq. 5 



II 
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At or above 3 Hz and SNR of 15 dB, batch processing produces range estimates of 
component motion and avoids the need for a phase tracking loop, since RMS error is less than 1 
radian. For a spin or precession component amplitude of .0075 m or less, no WBER spectrum 
detection is required. If the combined PDR and track filtering requires a phase tracking loop, 
options include processing at 3 Hz PRF and a SNR of 15 dB, with increase of data dwell, or 
increased chirp bandwidth, or both, or 6 Hz PRF and SNR of 12 dB for reduced error in 
estimating phase difference pulse-to-pulse. 

Although described in the context of a tracking radar system, it will be appreciated that 
the technique of the present invention could be used in other applications that require high- 
precision Doppler measurements (for example, weather radars). 

Other embodiments are within the scope of the following claims. 

What is claimed is: 
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